Dynamic recovery method and system for uavs and storage medium

ABSTRACT

The disclosure provides a dynamic recovery method and system for UAVs and storage medium. On the basis of obtaining the environmental information of UAVs, combined with different factors such as the following task quantity of each UAV, the distance from the current position of each UAV to the recovery point, current performance of each UAV, total cost function is constructed, and optimization is carried out through improved particle swarm optimization algorithm. Then on the premise of the lowest recovery cost, the priority of each UAV is obtained. For the ordered UAV, the route to the recovery point is determined. At the same time, it is necessary to consider conflict resolution of UAVs when encountering obstacles in the flight process, as well as the possible failure in line recovery. At this time, it is necessary to reorder some hovering UAVs to ensure flight safety of UAVs and robustness of the algorithm.

CROSS-REFERENCE TO RELATED APPLICATIONS

This application claims priority to Chinese Patent Application No. 202010196990.6, filed on Mar. 19, 2020, which is hereby incorporated by reference in its entirety.

TECHNICAL FIELD

The disclosure involves the technical field of UAV (Unmanned Aerial Vehicle), and relates to a dynamic recovery method and system of UAVs and storage medium.

BACKGROUND

UAV has many advantages such as low cost, low loss, high survival rate and so on. At present, it has been applied in many fields such as military reconnaissance, disaster monitoring, pesticide spraying, express transportation and so on. With the increasingly complex and changeable flight environment, the executive ability of a single UAV is difficult to meet the actual task requirements. Multi UAVs formation to complete the task together becomes the main mode in the future UAVs.

UAVs can control a larger range, achieve complementary functions and task sharing and so on. Therefore, UAVs can perform complex tasks, shorten the duration of performing task and improve the efficiency of attendance. It is of great significance for the UAVs formation to complete the given task cooperatively.

The task being completed cooperatively by UAVs cluster involves key technologies such as cluster control technology, cluster communication network design, task planning technology, cluster launch and recovery and so on. Among them, cluster launch and recovery technology are key steps in the task performed by the cluster. To improve the launch efficiency of UAVs means to form a large number of multi-wave UAVs flying rapidly. At the same time, as a simple rear support for UAVs flying, UAVs recovery technology can not only improve the quality of information obtained by the UAVs, but also reuse the recovered UAVs, greatly reducing UAVs cost.

Traditional UAV recovery technologies include parachute recovery technology, net recovery technology, line recovery technology and mechanical arm recovery technology and so on, and most of them are for a single UAV. With the development of UAV towards miniaturization and cluster, a large number of UAVs will be released in the first time when the UAVs are performing their tasks to establish the network in the air. The dynamic recovery of UAVs will be the simplest logistics support scheme for UAVs to perform tasks. How to dynamically recycle the entire UAVs is a problem to be solved.

SUMMARY

Based on the above problems, the disclosure provides a dynamic recovery method of UAVs, which ensures the robustness of the algorithm and realizes the dynamic recovery of UAVs by prioritizing the UAVs and flying them to the recovery point in order.

The dynamic recovery method of UAVs comprises the following steps:

step 1: in complex environment, UAVs obtain their own performance information and task information;

the task information includes: the current position of each UAV, the position of recovery point, the coordinates of obstacles within the detection radius of each UAV and the coordinate positions of other UAVs;

its performance information includes: the current flight time of UAV before recovery, the number of tasks to be performed by UAV and whether UAV is damaged.

Step 2: according to the performance information and task information of each UAV, the total cost function F of dynamic recovery for UAVs is established.

Firstly, for a UAV i, the total cost function ƒ_(i) is calculated according to its recovery order.

The calculation formula is as follows:

ƒ_(i)=ƒ_(1i)+ƒ_(2i)+ƒ_(3i)+ƒ_(4i)+ƒ_(5i)+ƒ_(6i)

ƒ_(1i) is the cost function related to the recovery order generated according to the number of tasks Q_(i) waiting to be completed of UAV i after recovery. The specific formula is as follows:

ƒ_(1i) =O _(i)·(Q _(i))^(α) ¹

wherein, O_(i) is a permutation of N UAVs; α₁ is the intensity parameter of the effect of the number of tasks waiting to be completed Q_(i) on the cost function;

ƒ_(2i) is the cost function related to the recovery order generated according to the time T_(i) ^(p) needed to complete other tasks of UAV i after recovery. The specific formula is as follows:

ƒ_(2i) =O _(i)·(T _(i) ^(p))^(α) ² ,

α₂ is the intensity parameter of the effect of the time needed to complete other tasks T_(i) ^(p) on the cost function.

ƒ_(3i) is the cost function related to the recovery order generated according to the distance L_(i) between the current position and the recovery point of UAV i. The specific formula is as follows:

${f_{3i} = \frac{\left( L_{i} \right)^{\alpha_{3}}}{O_{i}}},$

α₃ is the intensity parameter of the effect of the distance L_(i) on the cost function.

f_(4i) is the cost function related to the recovery order generated based on the number of obstacles C_(i) of UAV i that may be encountered in flight. The specific formula is as follows:

${f_{4i} = \frac{\left( C_{i} \right)^{\alpha_{4}}}{O_{i}}},$

α₄ is the intensity parameter of the effect of the number of obstacles C_(i) on the cost function.

ƒ_(5i) is the cost function related to the recovery order generated based on the current damage severity D_(i) of UAV i. The specific formula is as follows:

ƒ_(5i) =O _(i)·(D _(i))^(α) ⁵ ,

α₅ is the intensity parameter of the effect of the current damage severity D_(i) on the cost function.

ƒ_(6i) is the cost function related to the recovery order generated according to the current flight time T_(i) ^(r) of UAV i. The specific formula is as follows:

ƒ_(6i) =O _(i)·(T _(i) ^(r))^(α) ⁶

α₆ is the intensity parameter of the effect of the current flight time T_(i) ^(r) on the cost function.

Then, for N UAVs among the whole UAVs, the total cost function F generated by recovery order is as follows:

F=Σ _(i=1) ^(N)ƒ_(i),

After expansion, the cost function of the dynamic recovery for UAVs is obtained:

$F = {\sum\limits_{i = 1}^{N}\left( {{O_{i} \cdot \left( Q_{i} \right)^{\alpha_{1}}} + {O_{i} \cdot \left( T_{i}^{p} \right)^{\alpha_{2}}} + \frac{\left( L_{i} \right)^{\alpha_{3}}}{O_{i}} + \frac{\left( C_{i} \right)^{\alpha_{4}}}{O_{i}} + {O_{i} \cdot \left( D_{i} \right)^{\alpha_{5}}} + {O_{i} \cdot \left( T_{i}^{r} \right)^{\alpha_{6}}}} \right)}$

N is the number of all UAVs.

Step 3: The improved particle swarm optimization algorithm is used to optimize the total cost function F, and the optimal recovery order of N UAVs is obtained;

In particle swarm optimization algorithm, the speed and position updating formulas of the qth position of the sth particle in the (t+1) generation are as follows:

v _(sq)(t+1)=v _(sq)(t)+c ₁×rand( )×(pbest_(sq)(t)−x _(sq)(t))+c ₂×rand( )×(gbest_(sq)(t)−x _(sq)(t)),

x _(sq)(t+1)=x _(sq)(t)+v _(sq)(t+1).

Where s=1, 2, . . . , n, n is the number of particles. Each particle is a permutation of (1, 2, . . . , N), and each permutation corresponds to a value of the cost function F.

v_(sq)(t) is the velocity of the q th position of the sth particle in the t generation, c₁ and c₂ are two constants, and rand( ) is the uniform random number in [0,1], pbest_(sq)(t) is the qth position of the historical optimal solution of the sth particle in the t generation, x_(sq) (t) is the qth position of the sth particle in the t generation, gbest_(sq)(t) means that after the historical optimal solutions of the top three particles in the t generation are determined, the qth position corresponding to a solution is randomly selected from the three historical optimal solutions.

After the qth position of the sth particle is updated by the formula, the original position is replaced with the new position. All the initial particles are optimized by generation T, and the optimal solution is the optimal dynamic recovery order for UAVs.

Step 4: Each UAV is recovered according to the priority order, and whether the UAVs encounter conflict in the dynamic recovery process is judged. If yes, the UAV carries out conflict resolution and flies back in order. If not, the UAV flies back in order.

The conflict resolution calculation formula of UAV i is as follows:

${\overset{\rightarrow}{v_{I}}\left( {t + 1} \right)} = {\frac{{\overset{\rightarrow}{v_{\iota}}\left( {t + 1} \right)} + {w\overset{\rightarrow}{g_{\iota}}}}{{{\overset{\rightarrow}{v_{\iota}}\left( {t + 1} \right)} + {w\overset{\rightarrow}{g_{\iota}}}}}.}$

Where {right arrow over (v_(I))}(t+1) is the flight speed of UAV i in the next time step when encountering conflict, w is the influence weight of priority flight direction, {right arrow over (g_(l))} is the direction vector of priority flight of UAV i (from current position to the recovery point), {right arrow over (v_(l))}(t+1) is the adjustment speed of UAV i in the next time step when encountering conflict; the calculation formula is as follows:

${{\overset{\rightarrow}{v_{\iota}}\left( {t + 1} \right)} = {- {\sum_{j = 1}^{n_{D}}\frac{\overset{\rightarrow}{r_{ij}}(t)}{{\overset{\rightarrow}{r_{ij}}(t)}}}}};$

n_(D) is the number of obstacles and other UAVs within the detection radius R_(D) of UAV i at time t, and {right arrow over (r_(lJ))}(t) represents the relative position of UAV i, obstacles and other UAV j at time t.

Step 5: After some UAVs arrive at the specific recovery area, they wait in line and hover around to judge whether there are UAVs that fail to complete line recovery. If yes, we reorder some UAVs that failed in recovery. If not, we recycle in order until recovery is completed.

The cost formula F_(ƒ) of reordering is as follows:

F _(ƒ)=Σ_(i=1) ^(N) ^(ƒ) (O′ _(i)·(Q _(i))^(α) ¹ +O′ ₁·(T _(i) ^(p))^(α) ² +O′ _(i)·(D _(i))^(α) ⁵ +O′ _(i)·(T _(i) ^(r))^(α) ⁶ ).

O′_(i) is a permutation of N_(f) hovering UAVs.

The improved particle swarm optimization algorithm is used to solve the problem again, and the optimal solution is the optimal recovery order of UAVs hovering around.

The disclosure has the advantages of:

1. The disclosure provides a dynamic recovery method for UAVs, which can realize the efficient dynamic recovery of UAVs in complex environment, and further cooperate to complete the follow-up tasks, which provides a new solution for the dynamic recovery problem of UAVs.

2. The disclosure provides a dynamic recovery method for UAVs. Starting from improving the safety and efficiency of air traffic system, the dynamic recovery method for UAVs is studied to ensure the flight safety of UAV, reduce the flight cost and improve the use efficiency of UAVs.

3. The disclosure provides a dynamic recovery method for UAVs, which can deal with the dynamic recovery problem of a large number of UAVs, so that the UAV can quickly recycle and complete other tasks after completing the flight task. The algorithm has low complexity and high precision, which can effectively solve the dynamic recovery problem of UAVs in complex environment.

BRIEF DESCRIPTION OF DRAWINGS

FIG. 1 is a flow chart of a dynamic recovery method for UAVs according to the present application;

FIG. 2 is the schematic diagram of conflict resolution of UAV in the recovery process according to the present application.

DESCRIPTION OF EMBODIMENTS

In order to facilitate the understanding and implementation of the disclosure by ordinary technical personnel in the field, the disclosure is further described in detail and in depth in combination with the drawings.

The disclosure provides a dynamic recovery method for UAVs. On the basis of obtaining the environment information of UAVs, combining with factors such as the following task quantity of each UAV, the distance between each UAV and recovery point, the current performance of each UAV and the different environment and so on, the priority of UAVs is carried out by improved particle swarm optimization algorithm, which makes the total recovery cost of UAVs the lowest. After the recovery order is determined, the route of ordered UAVs to the recovery point should be determined. At the same time, if the UAVs encounter obstacles in the flight process, conflict resolution is needed to ensure the flight safety of UAVs. Finally, taking line recovery as an example, it is also necessary to consider the actual recovery failure of UAVs in line recovery. At this time, it is necessary to reorder some hovering UAVs to ensure the flight safety of UAVs and the robustness of the algorithm.

As shown in FIG. 1, the dynamic recovery method for UAVs includes the following steps:

step 1: UAVs obtain environment information and their own performance and task information.

In complex flight environment, UAV may face a variety of static or dynamic threats, such as trees, buildings, motor vehicles or other UAVs, so it is necessary to obtain environmental information.

Each UAV obtains its task information in the flight scene, including: the current position of each UAV, the position of UAV recovery point, so as to determine the flight direction of UAVs; it also includes the coordinates of obstacles within the detection radius of each UAV and the coordinates of other UAVs, so as to avoid obstacles in the flight process of each UAV;

UAV also needs to know its own performance and state information, such as the current flight time of each UAV before recovery, the number of tasks to be performed by each UAV after recovery, and whether the UAV is damaged or not.

Step 2: according to the performance information and task information of each UAV, the total cost function F of dynamic recovery for UAVs is established;

In the process of dynamic recovery of UAVs, we must consider various factors such as the different number of UAVs, the difference of task, the difference of flight environment, the current performance of UAVs and so on. The priority of dynamic recovery process of each UAV among UAVs can be ranked, and then the dynamic recovery can be carried out according to the priority. The optimization index is the result of priority ranking, that is, the recovery order of UAVs.

The disclosure considers the following factors for the recovery order of a single UAV, which are all cost functions related to the recovery order.

1. The task requirements of UAV.

Assuming that for UAV i, the number of tasks waiting to be completed after recovery is Q_(i) and the waiting time needed to complete the tasks is T_(i) ^(p), then the higher the Q_(i) and T_(i) ^(p), the more backward the recovery order of UAV i is, the greater the recovery order cost is. The cost of recovery order is positively correlated with the size of Q_(i) and T_(i) ^(p).

The recovery order cost of UAV i according to Q_(i) is defined as ƒ_(1i). The specific formula is as follows:

ƒ_(1i) =O _(i)·(Q _(i))^(α) ¹ .

Wherein, O_(i)(i=1, 2, . . . , N) is a permutation of N UAVs; α₁ is the intensity parameter of the effect of Q_(i) on the cost function. The reason why α₁ is placed in the power exponent is that the more tasks the UAV wait to complete, in addition to the time required to complete the tasks, additional time consumption such as recovery and maintenance should be considered.

Similarly, the recovery order cost of UAV i according to T_(i) ^(p) is defined as ƒ_(2i). The specific formula is as follows:

ƒ_(2i) =O _(i)·(T _(i) ^(p))^(α) ² .

Where α₂ is the intensity parameter of the effect of T_(i) ^(p) on the cost function.

2. The flight environment of UAV.

Assuming that the distance between the current position of UAV and the recovery point is L_(i) for UAV i, the environmental complexity of flying to the recovery point mainly considers the number of obstacles C_(i) that may be encountered. Different from the impact of task requirements on the cost, the higher the L_(i) and C_(i) are, the more backward the recovery order of UAV i is, the lower the recovery order cost of UAV i is. The cost of recovery order is inversely related to the size of L_(i) and C_(i). This is because the larger the L_(i) and C_(i), the more time it takes for UAV i to fly to the recovery point. If the recovery order is too forward, UAV i cannot fly to the recovery point in a short time. On the contrary, it will affect the recovery of other UAVs among the UAVs and increase the total recovery cost.

The recovery order cost of UAV i according to L_(i) is defined as ƒ_(3i). The specific formula is as follows:

$f_{3i} = \frac{\left( L_{i} \right)^{\alpha_{3}}}{O_{i}}$

Where α₃ is the intensity parameter of the effect of L_(i) on the cost function.

Similarly, the recovery order cost of UAV i according to C_(i) is defined as ƒ_(4i). The specific formula is as follows:

$f_{4i} = \frac{\left( C_{i} \right)^{\alpha_{4}}}{O_{i}}$

Where α₄ is the intensity parameter of the effect of C_(i) on the cost function.

3. The current performance of UAV.

Assuming that the current damage severity of UAV i is D_(i) and the current flight time of UAV i is T_(i) ^(r). Then, with the increase of D_(i) and T_(i) ^(r), the more backward the recovery order of UAV i is, the greater the recovery order cost will be. The cost of recovery order is positively correlated with the size of D_(i) and T_(i) ^(r). This is because the greater the damage UAV i suffers and the greater the current flight time (the shorter the remaining flight time), the more it needs to be recovered and maintained in time.

The recovery order cost of UAV i according to D_(i) is defined as ƒ_(5i). The specific formula is as follows:

ƒ_(5i) =O _(i)·(D _(i))^(α) ⁵

Where α₅ is the intensity parameter of the effect of D_(i) on the cost function.

Similarly, the recovery order cost of UAV i according to T_(i) ^(r) is defined as ƒ_(6i). The specific formula is as follows:

ƒ_(6i) =O _(i)·(T _(i) ^(r))^(α) ⁶

Where α₆ is the intensity parameter of the effect of T_(i) ^(r) on the cost function.

When the above factors are determined, the total cost of UAV i according to its recovery order O_(i) is defined as ƒ_(i), and the calculation formula of ƒ_(i) is as follows:

ƒ_(i)=ƒ_(1i)+ƒ_(2i)+ƒ_(3i)+ƒ_(4i)+ƒ_(5i)+ƒ_(6i)

Assuming that there are N UAVs in the whole UAV group, the total cost F generated by the whole UAV group according to its recovery order is as follows:

F=Σ _(i=1) ^(N)ƒ_(i).

After expansion, the formula is as follows:

$F = {\sum\limits_{i = 1}^{N}\left( {{O_{i} \cdot \left( Q_{i} \right)^{\alpha_{1}}} + {O_{i} \cdot \left( T_{i}^{p} \right)^{\alpha_{2}}} + \frac{\left( L_{i} \right)^{\alpha_{3}}}{O_{i}} + \frac{\left( C_{i} \right)^{\alpha_{4}}}{O_{i}} + {O_{i} \cdot \left( D_{i} \right)^{\alpha_{5}}} + {O_{i} \cdot \left( T_{i}^{r} \right)^{\alpha_{6}}}} \right)}$

Accordingly, the overall recovery cost corresponding to the priority sequence of the whole UAV group is obtained. The lower the cost, the higher the profit. And the next goal is to optimize the value of the independent variables of the objective function, that is, the final recovery order of UAVs, so that the overall recovery cost is low.

Step 3: The improved particle swarm optimization algorithm is used to optimize the total cost function F, and the optimal recovery order of UAVs is obtained;

Initialize n recovery orders, each recovery order is a permutation of (1, 2, . . . , N) and corresponds to a value of the cost function F. The goal is to iteratively optimize the initial recovery order through particle swarm optimization algorithm, and each recovery order is constantly updated. After a certain number of iterations, an optimal recovery order is finally obtained, that is to make the value of cost function F the lowest.

In particle swarm optimization algorithm, the speed and position updating formula of the qth position of the sth particle (that is the sth recovery order) in the (t+1) generation are as follows:

v _(sq)(t+1)=v _(sq)(t)+c ₁×rand( )×(pbest_(sq)(t)−x _(sq)(t))+c ₂×rand( )×(gbest_(sq)(t)−x _(sq)(t)).

x _(sq)(t+1)=x _(sq)(t)+v _(sq)(t+1).

Where s=1, 2, . . . , n, v_(sq)(t) is the velocity of the qth position of the sth particle (that is the sth recovery order) in the t generation, x_(sq)(t) is the qth position of the sth particle (that is the sth recovery order) in the t generation. The recovery order of each particle is evaluated by the cost function. c₁ and c₂ are two constants, and rand( ) is the uniform random number in [0,1], pbest_(sq)(t) is the qth position of the historical optimal solution of the sth particle in the t generation, gbest_(sq)(t) means that after the historical optimal solutions of the top three particles in the t generation are determined, the qth position corresponding to a solution is randomly selected from the three historical optimal solutions. The reason why we don't choose the best solution, but choose the best solution in the first three randomly is to ensure the randomness of the algorithm.

It is worth noting that after the qth position of the sth particle (that is the sth recovery order) is updated by formula, it is necessary to make improvement based on the background of the problem which is different from the original PSO algorithm directly updating the next position. Because each solution is a permutation of (1, 2, . . . , N), after updating a certain position, the original position must be replaced with the new position. This is to ensure that the updated solution is also a permutation of (1, 2, . . . , N) and the algorithm can continue. After all the initial particles are optimized by iteration of generation T, the final optimal solution is the optimal recovery order of UAVs.

Step 4: Each UAV is recovered according to the priority order, and whether the UAVs encounters conflict in the dynamic recovery process is judged. If yes, the UAV carries out conflict resolution and flies back in order. If not, the UAV flies back in order.

As shown in FIG. 2, the UAV has a detection radius R_(D). Taking UAV i as an example, it is assumed that the number of obstacles and other UAVs within the detection radius of UAV i at time t is n_(D), then the flight speed of UAV i in the next time step is adjusted according to the following formula:

${\overset{\rightarrow}{v_{\iota}}\left( {t + 1} \right)} = {- {\sum_{j = 1}^{n_{D}}{\frac{\overset{\rightarrow}{r_{ij}}(t)}{{\overset{\rightarrow}{r_{ij}}(t)}}.}}}$

Where {right arrow over (r_(lJ))}(t) is the relative position of UAV i, obstacles and other UAV j at time t, i.e. UAV i flies away from risk, which is of great significance for the flight safety of UAV.

At the same time, UAV i also has a priority flight direction vector {right arrow over (g_(l))}. the direction of vector {right arrow over (g_(l))} is the direction that the current position of UAV i points to the recovery point.

Considering the above two factors, when the UAV i encounters a conflict, the flight speed of UAV i in the next time step is:

${{\overset{\rightarrow}{v}}_{I}\left( {t + 1} \right)} = {\frac{{\overset{\rightarrow}{v_{\iota}}\left( {t + 1} \right)} + {w\;\overset{\rightarrow}{g_{\iota}}}}{{{\overset{\rightarrow}{v_{\iota}}\left( {t + 1} \right)} + {w\;\overset{\rightarrow}{g_{\iota}}}}}.}$

Where w is the weight of the priority direction. In this way, the UAV i can not only fly to the target direction, but also achieve conflict resolution.

Step 5: After some UAVs arrive at the specific recovery area, they wait in line and hover around to judge whether there are UAVs that fail to complete line recovery. If yes, we reorder some UAVs that failed in recovery. If not, we recycle in order until recovery is completed.

The UAV may encounter the situation of actual recovery failure in line recovery. At this time, it is necessary to reorder some of the UAVs hovering and waiting (assuming its number is N_(ƒ)) among the UAVs to ensure the robustness of the algorithm.

At this time, because the UAV is near the recovered ship, the UAVs flight environment is not considered, only the UAVs task requirements and current performance are considered, and the UAVs that failed to be recovered and some of the UAVs hovering and waiting are reordered. The cost formula F_(ƒ) is defined as follows:

F _(ƒ)=Σ_(i=1) ^(N) ^(ƒ) (O′ _(i)·(Q _(i))^(α) ¹ +O′ ₁·(T _(i) ^(p))^(α) ² +O′ _(i)·(D _(i))^(α) ⁵ +O′ _(i)·(T _(i) ^(r))^(α) ⁶ ).

Where O′(i=1, 2, . . . , N_(ƒ)) is a permutation of (1, 2, . . . , N_(ƒ)).

The improved particle swarm optimization algorithm is used to solve the problem again, and the optimal solution is the optimal recovery order of UAVs hovering and waiting.

The disclosure can realize the dynamic recovery of UAVs under complex conditions, realize the efficient flight and dynamic recovery of UAVs, and complete the assigned tasks cooperatively, which has positive significance.

A storage medium comprising instructions that when executed cause a processor of a computing device to:

in complex environment, obtain their own performance information and task information by UAVs;

wherein the task information comprises: the current position of each UAV, the position of recovery point, the coordinates of obstacles within the detection radius of each UAV and the coordinate positions of other UAVs; its performance information comprises: the current flight time of each UAV before recovery, the number of tasks to be performed by each UAV and whether the UAV is damaged,

according to the performance information and task information of each UAV, establish the total cost function F of dynamic recovery for UAVs, and the total cost function is calculated as follows:

$F = {\sum\limits_{i = 1}^{N}\left( {{O_{i} \cdot \left( Q_{i} \right)^{\alpha_{1}}} + {O_{i} \cdot \left( T_{i}^{p} \right)^{\alpha_{2}}} + \frac{\left( L_{i} \right)^{\alpha_{3}}}{O_{i}} + \frac{\left( C_{i} \right)^{\alpha_{4}}}{O_{i}} + {O_{i} \cdot \left( D_{i} \right)^{\alpha_{5}}} + {O_{i} \cdot \left( T_{i}^{r} \right)^{\alpha_{6}}}} \right)}$

wherein, O_(i) is a permutation of N UAVs; Q_(i) is the number of tasks waiting to be completed after UAV i is recovered; α₁ is the intensity parameter of the effect of the number of tasks waiting to be completed Q_(i) on the cost function; T_(i) ^(p) is the waiting time needed to complete other tasks after UAV i is recovered; α₂ is the intensity parameter of the effect of the waiting time needed to complete other tasks T_(i) ^(p) on the cost function; L_(i) is the distance from the current position of UAV i to the recovery point; α₃ is the intensity parameter of the effect of the distance L_(i) on the cost function; C_(i) is the number of obstacles that UAV i may encounter in flight; α₄ is the intensity parameter of the effect of the number of obstacles C_(i) on the cost function; D_(i) is the current damage severity of UAV i; α₅ is the intensity parameter of the effect of the current damage severity D_(i) on the cost function; T_(i) ^(r) is the current flight time of UAV i; α₆ is the intensity parameter of the effect of the current flight time T_(i) ^(r) on the cost function;

by using the improved particle swarm optimization algorithm, optimize the total cost function F, and obtain the optimal recovery order of UAVs;

in particle swarm optimization algorithm, the speed and position updating formula of the qth position of the sth particle in the (t+1) generation are as follows:

v _(sq)(t+1)=v _(sq)(t)+c ₁×rand( )×(pbest_(sq)(t)−x _(sq)(t))+c ₂×rand( )×(gbest_(sq)(t)−x _(sq)(t));

x _(sq)(t+1)=x _(sq)(t)+v _(sq)(t+1);

wherein s=1, 2, . . . , n is the number of particles, each particle is a permutation of (1, 2, . . . , N), and each permutation corresponds to a value of the cost function F;

v_(sq)(t) is the velocity of the qth position of the sth particle in the t generation, c_(t) and c₂ are two constants, and rand( ) is the uniform random number in [0,1], pbest_(sq)(t) is the qth position of the historical optimal solution of the sth particle in the t generation, x_(sq) (t) is the qth position of the sth particle in the t generation, gbest_(sq)(t) means that after the historical optimal solutions of the top three particles in the t generation are determined, the qth position corresponding to a solution is randomly selected from the three historical optimal solutions;

recover each UAV according to the priority order, and judge whether the UAVs encounter conflict in the dynamic recovery process, if yes, the UAV carries out conflict resolution and flies back in order, if not, the UAV flies back in order,

the conflict resolution calculation formula of UAV i is as follows:

${{{\overset{\rightarrow}{v}}_{I}\left( {t + 1} \right)} = \frac{{\overset{\rightarrow}{v_{\iota}}\left( {t + 1} \right)} + {w\;\overset{\rightarrow}{g_{\iota}}}}{{{\overset{\rightarrow}{v_{\iota}}\left( {t + 1} \right)} + {w\;\overset{\rightarrow}{g_{\iota}}}}}},$

where {right arrow over (v_(I))}(t+1) is the flight speed of UAV i in the next time step when encountering conflict, w is the influence weight of priority flight direction, {right arrow over (g_(l))} is the direction vector of priority flight of UAV i (from current position to the recovery point); {right arrow over (v_(l))}(t+1) is the adjustment speed of UAV i in the next time step when encountering conflict;

judge whether there are UAVs that fail to complete line recovery, after some UAVs arriving at the specific recovery area, if yes, we reorder some UAVs that failed in recovery, if not, we recycle in order until recovery is completed.

A dynamic recovery system for UAVs, comprising: the UAVs, a first processor, a second processor, a third processor, and a fourth processor, wherein

in complex environment, UAVs obtain their own performance information and task information;

a first processor, according to the performance information and task information of each UAV, establishes the total cost function F of dynamic recovery for UAVs, and calculates the total cost function as follows:

$F = {\sum\limits_{i = 1}^{N}\left( {{O_{i} \cdot \left( Q_{i} \right)^{\alpha_{1}}} + {O_{i} \cdot \left( T_{i}^{p} \right)^{\alpha_{2}}} + \frac{\left( L_{i} \right)^{\alpha_{3}}}{O_{i}} + \frac{\left( C_{i} \right)^{\alpha_{4}}}{O_{i}} + {O_{i} \cdot \left( D_{i} \right)^{\alpha_{5}}} + {O_{i} \cdot \left( T_{i}^{r} \right)^{\alpha_{6}}}} \right)}$

wherein, O_(i) is a permutation of N UAVs; Q_(i) is the number of tasks waiting to be completed after UAV i is recovered; α₁ is the intensity parameter of the effect of the number of tasks waiting to be completed Q_(i) on the cost function; T_(i) ^(p) is the waiting time needed to complete other tasks after UAV i is recovered; α₂ is the intensity parameter of the effect of the waiting time needed to complete other tasks T_(i) ^(p) on the cost function; L_(i) is the distance from the current position of UAV i to the recovery point; α₃ is the intensity parameter of the effect of the distance L_(i) on the cost function; C_(i) is the number of obstacles that UAV i may encounter in flight; α₄ is the intensity parameter of the effect of the number of obstacles C_(i) on the cost function; D_(i) is the current damage severity of UAV i; α₅ is the intensity parameter of the effect of the current damage severity D_(i) on the cost function; T_(i) ^(r) is the current flight time of UAV i; α₆ is the intensity parameter of the effect of the current flight time T_(i) ^(r) on the cost function;

the second processor uses the improved particle swarm optimization algorithm to optimize the total cost function F, and obtains the optimal recovery order of UAVs;

in particle swarm optimization algorithm, the speed and position updating formula of the qth position of the sth particle in the (t+1) generation are as follows:

v _(sq)(t+1)=v _(sq)(t)+c ₁×rand( )×(pbest_(sq)(t)−x _(sq)(t))+c ₂×rand( )×(gbest_(sq)(t)−x _(sq)(t)),

x _(sq)(t+1)=x _(sq)(t)+v _(sq)(t+1);

wherein s=1, 2, . . . , n, n is the number of particles, each particle is a permutation of (1, 2, . . . , N), and each permutation corresponds to a value of the cost function F;

v_(sq)(t) is the velocity of the qth position of the sth particle in the t generation, c₁ and c₂ are two constants, and rand( ) is the uniform random number in [0,1], pbest_(sq)(t) is the qth position of the historical optimal solution of the sth particle in the t generation, x_(sq) (t) is the qth position of the sth particle in the t generation, gbest_(sq)(t) means that after the historical optimal solutions of the top three particles in the t generation are determined, the qth position corresponding to a solution is randomly selected from the three historical optimal solutions;

the third processor recovers each UAV according to the priority order, and judges whether the UAVs encounter conflict in the dynamic recovery process, if yes, the UAV carries out conflict resolution and flies back in order, if not, the UAV flies back in order,

the conflict resolution calculation formula of UAV i is as follows:

${{{\overset{\rightarrow}{v}}_{I}\left( {t + 1} \right)} = \frac{{\overset{\rightarrow}{v_{\iota}}\left( {t + 1} \right)} + {w\;\overset{\rightarrow}{g_{\iota}}}}{{{\overset{\rightarrow}{v_{\iota}}\left( {t + 1} \right)} + {w\;\overset{\rightarrow}{g_{\iota}}}}}},$

where {right arrow over (v_(I))}(t+1) is the flight speed of UAV i in the next time step when encountering conflict, w is the influence weight of priority flight direction, {right arrow over (g_(l))} is the direction vector of priority flight of UAV i (from current position to the recovery point); {right arrow over (v_(l))}(t+1) is the adjustment speed of UAV i in the next time step when encountering conflict;

The formula for adjustment speed is:

${{\overset{\rightarrow}{v_{\iota}}\left( {t + 1} \right)} = {- {\sum\limits_{j = 1}^{n_{D}}\frac{\overset{\rightarrow}{r_{\iota\; J}}(t)}{{\overset{\rightarrow}{r_{\iota\; J}}(t)}}}}},$

n_(D) is the number of obstacles and other UAVs within the detection radius R_(D) of UAV i at time t, and {right arrow over (r_(lJ))}(t) represents the relative position of UAV i, obstacles and other UAV j at time t.

the fourth processor judges whether there are UAVs that fail to complete line recovery when they wait in line and hover around, after some UAVs arrive at the specific recovery area, if yes, we reorder some UAVs that failed in recovery, if not, we recycle in order until recovery is completed.

The cost formula F_(ƒ) of reordering is as follows:

$F_{f} = {\sum\limits_{i = 1}^{N_{f}}\left( {{O_{i}^{\prime} \cdot \left( Q_{i} \right)^{\alpha_{1}}} + {O_{i}^{\prime} \cdot \left( T_{i}^{p} \right)^{\alpha_{2}}} + {O_{i}^{\prime} \cdot \left( D_{i} \right)^{\alpha_{5}}} + {O_{i}^{\prime} \cdot \left( T_{i}^{r} \right)^{\alpha_{6}}}} \right)}$

O′_(i) is a permutation of N_(ƒ) hovering UAVs;

the improved particle swarm optimization algorithm is used to solve the problem again, and the optimal solution is the optimal recovery order of UAVs hovering around.

Each of the first processor, the second processor, the third processor and the fourth processor is independent processor, or all of them are integrated in a single processor.

The foregoing descriptions of specific exemplary embodiments of the present application have been presented for purposes of illustration and description. They are not intended to be exhaustive or to limit the application to the precise forms disclosed, and obviously many modifications and variations are possible in light of the above teachings. The exemplary embodiments were chosen and described in order to explain certain principles of the application and their practical application, to thereby enable others skilled in the art to make and utilize various exemplary embodiments of the present application, as well as various alternatives and modifications thereof. It is intended that the scope of the application be defined by the claims appended hereto and their equivalents. 

What is claimed is:
 1. A dynamic recovery method for UAVs, comprising the following steps: step 1: in complex environment, UAVs obtain their own performance information and task information; step 2: according to the performance information and task information of each UAV, the total cost function F of dynamic recovery for UAVs is established, and the total cost function is calculated as follows: $F = {\sum\limits_{i = 1}^{N}\left( {{O_{i} \cdot \left( Q_{i} \right)^{\alpha_{1}}} + {O_{i} \cdot \left( T_{i}^{p} \right)^{\alpha_{2}}} + \frac{\left( L_{i} \right)^{\alpha_{3}}}{O_{i}} + \frac{\left( C_{i} \right)^{\alpha_{4}}}{O_{i}} + {O_{i} \cdot \left( D_{i} \right)^{\alpha_{5}}} + {O_{i} \cdot \left( T_{i}^{r} \right)^{\alpha_{6}}}} \right)}$ wherein, O_(i) is a permutation of N UAVs; Q_(i) is the number of tasks waiting to be completed after UAV i is recovered; α₁ is the intensity parameter of the effect of the number of tasks waiting to be completed Q_(i) on the cost function; T_(i) ^(p) is the waiting time needed to complete other tasks after UAV i is recovered; α₂ is the intensity parameter of the effect of the waiting time needed to complete other tasks T_(i) ^(p) on the cost function; L_(i) is the distance from the current position of UAV i to the recovery point; α₃ is the intensity parameter of the effect of the distance L_(i) on the cost function; C_(i) is the number of obstacles that UAV i may encounter in flight; α₄ is the intensity parameter of the effect of the number of obstacles C_(i) on the cost function; D_(i) is the current damage severity of UAV i; α₅ is the intensity parameter of the effect of the current damage severity D_(i) on the cost function; T_(i) ^(r) is the current flight time of UAV i; α₆ is the intensity parameter of the effect of the current flight time T_(i) ^(r) on the cost function; step 3: the improved particle swarm optimization algorithm is used to optimize the total cost function F, and the optimal recovery order of UAVs is obtained; in particle swarm optimization algorithm, the speed and position updating formula of the qth position of the sth particle in the (t+1) generation are as follows: v _(sq)(t+1)=v _(sq)(t)+c ₁×rand( )×(pbest_(sq)(t)−x _(sq)(t))+c ₂×rand( )×(gbest_(sq)(t)−x _(sq)(t)), x _(sq)(t+1)=x _(sq)(t)+v _(sq)(t+1); wherein s=1, 2, . . . , n, n is the number of particles, each particle is a permutation of (1, 2, . . . , N), and each permutation corresponds to a value of the cost function F; v_(sq)(t) is the velocity of the qth position of the sth particle in the t generation, c₁ and c₂ are two constants, and rand( ) is the uniform random number in [0,1], pbest_(sq)(t) is the qth position of the historical optimal solution of the sth particle in the t generation, x_(sq)(t) is the qth position of the sth particle in the t generation, gbest_(sq)(t) means that after the historical optimal solutions of the top three particles in the t generation are determined, the q th position corresponding to a solution is randomly selected from the three historical optimal solutions; step 4: each UAV is recovered according to the priority order, and whether the UAVs encounter conflict in the dynamic recovery process is judged, if yes, the UAV carries out conflict resolution and flies back in order, if not, the UAV flies back in order, the conflict resolution calculation formula of UAV i is as follows: ${{\overset{\rightarrow}{v_{I}}\left( {t + 1} \right)} = \frac{{\overset{\rightarrow}{v_{\iota}}\left( {t + 1} \right)} + {w\;\overset{\rightarrow}{g_{\iota}}}}{{{\overset{\rightarrow}{v_{\iota}}\left( {t + 1} \right)} + {w\;\overset{\rightarrow}{g_{\iota}}}}}},$ where {right arrow over (v_(I))}(t+1) is the flight speed of UAV i in the next time step when encountering conflict, w is the influence weight of priority flight direction, {right arrow over (g_(l))} is the direction vector of priority flight of UAV i (from current position to the recovery point); {right arrow over (v_(l))}(t+1) is the adjustment speed of UAV i in the next time step when encountering conflict; step 5: after some UAVs arrive at the specific recovery area, they wait in line and hover around to judge whether there are UAVs that fail to complete line recovery, if yes, we reorder some UAVs that failed in recovery, if not, we recycle in order until recovery is completed.
 2. The dynamic recovery method for UAVs according to claim 1, wherein the task information described in step 1 comprises: the current position of each UAV, the position of recovery point, the coordinates of obstacles within the detection radius of each UAV and the coordinate positions of other UAVs; its performance information comprises: the current flight time of each UAV before recovery, the number of tasks to be performed by each UAV and whether the UAV is damaged.
 3. The dynamic recovery method for UAVs according to claim 1, the formula for adjustment speed in step 4 is: ${{\overset{\rightarrow}{v_{\iota}}\left( {t + 1} \right)} = {- {\sum\limits_{j = 1}^{n_{D}}\frac{\overset{\rightarrow}{r_{\iota\; J}}(t)}{{\overset{\rightarrow}{r_{\iota\; J}}(t)}}}}},$ n_(D) is the number of obstacles and other UAVs within the detection radius R_(D) of UAV i at time t, and {right arrow over (r_(lJ))}(t) represents the relative position of UAV i, obstacles and other UAV j at time t.
 4. The dynamic recovery method for UAVs according to claim 1, wherein the cost formula F_(ƒ) of reordering in step 5 is as follows: F _(ƒ)=Σ_(i=1) ^(N) ^(ƒ) (O′ _(i)·(Q _(i))^(α) ¹ +O′ ₁·(T _(i) ^(p))^(α) ² +O′ _(i)·(D _(i))^(α) ⁵ +O′ _(i)·(T _(i) ^(r))^(α) ⁶ ), O′_(i) is a permutation of N_(f) hovering UAVs; the improved particle swarm optimization algorithm is used to solve the problem again, and the optimal solution is the optimal recovery order of UAVs hovering around.
 5. A storage medium comprising instructions that when executed cause a processor of a computing device to: in complex environment, obtain their own performance information and task information by UAVs; according to the performance information and task information of each UAV, establish the total cost function F of dynamic recovery for UAVs, and the total cost function is calculated as follows: $F = {\sum\limits_{i = 1}^{N}\left( {{O_{i} \cdot \left( Q_{i} \right)^{\alpha_{1}}} + {O_{i} \cdot \left( T_{i}^{p} \right)^{\alpha_{2}}} + \frac{\left( L_{i} \right)^{\alpha_{3}}}{O_{i}} + \frac{\left( C_{i} \right)^{\alpha_{4}}}{O_{i}} + {O_{i} \cdot \left( D_{i} \right)^{\alpha_{5}}} + {O_{i} \cdot \left( T_{i}^{r} \right)^{\alpha_{6}}}} \right)}$ wherein, O_(i) is a permutation of N UAVs; Q_(i) is the number of tasks waiting to be completed after UAV i is recovered; α₁ is the intensity parameter of the effect of the number of tasks waiting to be completed Q_(i) on the cost function; T_(i) ^(p) is the waiting time needed to complete other tasks after UAV i is recovered; α₂ is the intensity parameter of the effect of the waiting time needed to complete other tasks T_(i) ^(p) on the cost function; L_(i) is the distance from the current position of UAV i to the recovery point; α₃ is the intensity parameter of the effect of the distance L_(i) on the cost function; C_(i) is the number of obstacles that UAV i may encounter in flight; α₄ is the intensity parameter of the effect of the number of obstacles C_(i) on the cost function; D_(i) is the current damage severity of UAV i; α₅ is the intensity parameter of the effect of the current damage severity D_(i) on the cost function; T_(i) ^(r) is the current flight time of UAV i; α₆ is the intensity parameter of the effect of the current flight time T_(i) ^(r) on the cost function; by using the improved particle swarm optimization algorithm, optimize the total cost function F, and obtain the optimal recovery order of UAVs; in particle swarm optimization algorithm, the speed and position updating formula of the qth position of the sth particle in the (t+1) generation are as follows: v _(sq)(t+1)=v _(sq)(t)+c ₁×rand( )×(pbest_(sq)(t)−x _(sq)(t))+c ₂×rand( )×(gbest_(sq)(t)−x _(sq)(t)), x _(sq)(t+1)=x _(sq)(t)+v _(sq)(t+1); wherein s=1, 2, . . . , n, n is the number of particles, each particle is a permutation of (1, 2, . . . , N), and each permutation corresponds to a value of the cost function F; v_(sq)(t) is the velocity of the qth position of the sth particle in the t generation, c₁ and c₂ are two constants, and rand( ) is the uniform random number in [0,1], pbest_(sq)(t) is the qth position of the historical optimal solution of the sth particle in the t generation, x_(sq)(t) is the qth position of the sth particle in the t generation, gbest_(sq)(t) means that after the historical optimal solutions of the top three particles in the t generation are determined, the qth position corresponding to a solution is randomly selected from the three historical optimal solutions; recover each UAV according to the priority order, and judge whether the UAVs encounter conflict in the dynamic recovery process, if yes, the UAV carries out conflict resolution and flies back in order, if not, the UAV flies back in order, the conflict resolution calculation formula of UAV i is as follows: ${{\overset{\rightarrow}{v_{I}}\left( {t + 1} \right)} = \frac{{\overset{\rightarrow}{v_{\iota}}\left( {t + 1} \right)} + {w\;\overset{\rightarrow}{g_{\iota}}}}{{{\overset{\rightarrow}{v_{\iota}}\left( {t + 1} \right)} + {w\;\overset{\rightarrow}{g_{\iota}}}}}},$ where {right arrow over (v_(I))}(t+1) is the flight speed of UAV i in the next time step when encountering conflict, w is the influence weight of priority flight direction, {right arrow over (g_(l))} is the direction vector of priority flight of UAV i (from current position to the recovery point); {right arrow over (v_(l))}(t+1) is the adjustment speed of UAV i in the next time step when encountering conflict; judge whether there are UAVs that fail to complete line recovery, after some UAVs arriving at the specific recovery area, if yes, we reorder some UAVs that failed in recovery, if not, we recycle in order until recovery is completed.
 6. A dynamic recovery system for UAVs, comprising: the UAVs, a first processor, a second processor, a third processor, and a fourth processor, wherein in complex environment, UAVs obtain their own performance information and task information; a first processor, according to the performance information and task information of each UAV, establishes the total cost function F of dynamic recovery for UAVs, and calculates the total cost function as follows: $F = {\sum\limits_{i = 1}^{N}\left( {{O_{i} \cdot \left( Q_{i} \right)^{\alpha_{1}}} + {O_{i} \cdot \left( T_{i}^{p} \right)^{\alpha_{2}}} + \frac{\left( L_{i} \right)^{\alpha_{3}}}{O_{i}} + \frac{\left( C_{i} \right)^{\alpha_{4}}}{O_{i}} + {O_{i} \cdot \left( D_{i} \right)^{\alpha_{5}}} + {O_{i} \cdot \left( T_{i}^{r} \right)^{\alpha_{6}}}} \right)}$ wherein, O_(i) is a permutation of N UAVs; Q_(i) is the number of tasks waiting to be completed after UAV i is recovered; α₁ is the intensity parameter of the effect of the number of tasks waiting to be completed Q_(i) on the cost function; T_(i) ^(p) is the waiting time needed to complete other tasks after UAV i is recovered; α₂ is the intensity parameter of the effect of the waiting time needed to complete other tasks T_(i) ^(p) on the cost function; L_(i) is the distance from the current position of UAV i to the recovery point; α₃ is the intensity parameter of the effect of the distance L_(i) on the cost function; C_(i) is the number of obstacles that UAV i may encounter in flight; α₄ is the intensity parameter of the effect of the number of obstacles C_(i) on the cost function; D_(i) is the current damage severity of UAV i; α₅ is the intensity parameter of the effect of the current damage severity D_(i) on the cost function; T_(i) ^(r) is the current flight time of UAV i; α₆ is the intensity parameter of the effect of the current flight time T_(i) ^(r) on the cost function; the second processor uses the improved particle swarm optimization algorithm to optimize the total cost function F, and obtains the optimal recovery order of UAVs; in particle swarm optimization algorithm, the speed and position updating formula of the qth position of the sth particle in the (t+1) generation are as follows: v _(sq)(t+1)=v _(sq)(t)+c ₁×rand( )×(pbest_(sq)(t)−x _(sq)(t))+c ₂×rand( )×(gbest_(sq)(t)−x _(sq)(t)); x _(sq)(t+1)=x _(sq)(t)+v _(sq)(t+1); wherein s=1, 2, . . . , n, n is the number of particles, each particle is a permutation of (1, 2, . . . , N), and each permutation corresponds to a value of the cost function F; v_(sq)(t) is the velocity of the qth position of the sth particle in the t generation, c₁ and c₂ are two constants, and rand( ) is the uniform random number in [0,1], pbest_(sq)(t) is the qth position of the historical optimal solution of the sth particle in the t generation, x_(sq)(t) is the qth position of the sth particle in the t generation, gbest_(sq)(t) means that after the historical optimal solutions of the top three particles in the t generation are determined, the q th position corresponding to a solution is randomly selected from the three historical optimal solutions; the third processor recovers each UAV according to the priority order, and judges whether the UAVs encounter conflict in the dynamic recovery process, if yes, the UAV carries out conflict resolution and flies back in order, if not, the UAV flies back in order, the conflict resolution calculation formula of UAV i is as follows: ${{\overset{\rightarrow}{v_{I}}\left( {t + 1} \right)} = \frac{{\overset{\rightarrow}{v_{\iota}}\left( {t + 1} \right)} + {w\;\overset{\rightarrow}{g_{\iota}}}}{{{\overset{\rightarrow}{v_{\iota}}\left( {t + 1} \right)} + {w\;\overset{\rightarrow}{g_{\iota}}}}}},$ where {right arrow over (v_(I))}(t+1) is the flight speed of UAV i in the next time step when encountering conflict, w is the influence weight of priority flight direction, {right arrow over (g_(l) )} is the direction vector of priority flight of UAV i (from current position to the recovery point); {right arrow over (v_(l))}(t+1) is the adjustment speed of UAV i in the next time step when encountering conflict; the fourth processor judges whether there are UAVs that fail to complete line recovery when they wait in line and hover around after some UAVs arrive at the specific recovery area, if yes, we reorder some UAVs that failed in recovery, if not, we recycle in order until recovery is completed.
 7. The dynamic recovery system for UAVs according to claim 6, wherein the task information comprises: the current position of each UAV, the position of recovery point, the coordinates of obstacles within the detection radius of each UAV and the coordinate positions of other UAVs; its performance information comprises: the current flight time of each UAV before recovery, the number of tasks to be performed by each UAV and whether the UAV is damaged.
 8. The dynamic recovery system for UAVs according to claim 6, the formula for adjustment speed is: ${{\overset{\rightarrow}{v_{\iota}}\left( {t + 1} \right)} = {- {\sum\limits_{j = 1}^{n_{D}}\frac{\overset{\rightarrow}{r_{\iota\; J}}(t)}{{\overset{\rightarrow}{r_{\iota\; J}}(t)}}}}},$ n_(D) is the number of obstacles and other UAVs within the detection radius R_(D) of UAV i at time t, and {right arrow over (r_(lJ))}(t) represents the relative position of UAV i, obstacles and other UAV j at time t.
 9. The dynamic recovery system for UAVs according to claim 6, wherein the cost formula F_(ƒ) of reordering is as follows: F _(ƒ)=Σ_(i=1) ^(N) ^(ƒ) (O′ _(i)·(Q _(i))^(α) ¹ +O′ ₁·(T _(i) ^(p))^(α) ² +O′ _(i)·(D _(i))^(α) ⁵ +O′ _(i)·(T _(i) ^(r))^(α) ⁶ ), O′_(i) is a permutation of N_(ƒ) hovering UAVs; the improved particle swarm optimization algorithm is used to solve the problem again, and the optimal solution is the optimal recovery order of UAVs hovering around. 